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Abstract 

' We analyze the cohective modes of high-temperature QCD in the case when there is an anisotropy 

' in the momentum-space distribution function for the gluons. We perform a tensor decomposition of 

the gluon self-energy and solve the dispersion relations for both stable and unstable modes. Results 
are presented for a class of anisotropic distribution functions which can be obtained by stretching 
or squeezing an isotropic distribution function along one direction in momentum space. We find 
that there are three stable modes and either one or two unstable modes depending on whether 
the distribution function is stretched or squeezed. The presence of unstable modes which have 
exponential growth can lead to a more rapid thermalization and isotropization of the soft modes 
, in a quark gluon plasma and therefore may play an important role in the dynamical evolution of 

^ I a quark-gluon plasma. 
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I. INTRODUCTION 



In the ongoing ultra-relativistic heavy collision experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) and the upcoming ones at the Large Hadron Collider (LHC) the behavior of nuclear 
matter under extreme conditions will be studied. The hope of these experiments is to create 
temperatures which are high enough for nuclear matter to undergo a phase transition to a quark 
gluon plasma (QGP). The quark gluon plasma, if generated, is expected to expand, cool, and then 
hadronize in the final stage of its evolution. In this context, an outstanding question faced by 
experimentalists and theorists is whether or not the system will "thermalize" fast enough to allow 
a thermodynamic description of the system during the central part of its evolution. 

In this paper we study the role of the collective modes of finite-temperature QCD in the thermal- 
ization, particularly the isotropization, of a finite-temperature QGP with anisotropic momentum- 
space distribution functions. This question has been addressed in previous papers in which the ex- 
istence of instabilities of a QGP were studied. In Refs. 0,1^ Mrowczyhski discussed the existence 
of instabilities to chromomagnetic fluctuations with a particular orientation of the chromoelectric 
field and wave- vector. In those papers Mrowczyhski showed that there existed an instability which 
was the equivalent of the Weibel or filamentation instability in electrodynamics Q . Weibel showed 
in his original paper that, within electrodynamics, unstable transverse modes exist in plasmas with 
anisotropic momentum distributions and he also derived their rate of growth in linear response 
theory. These types of instabilities are potentially very important to QGP evolution at RHIC or 
LHC due to the large amount of momentum-space anisotropy in the gluon distribution functions 
at r ~ 1 fm/c. 

Mrowczyhski and Randrup have recently performed phenomenological estimates of the growth 
rate of the instabilities for two types of anisotropic distribution functions J^. They found that the 
degree of amplification of the Weibel instability is not expected to dominate the dynamics of a QGP 
instead being comparable to the contribution from elastic Boltzmann collisions. However, they did 
point out that since a large number of the unstable modes could be excited then it is possible that 
their combined effect on the overall dynamics could be significant. In this paper we perform a 
detailed study of the hard-thermal-loop resummed gluon self-energy including a complete tensor 
decomposition of the self-energy and identification of all stable and unstable collective modes. 

In Sec.niwe set up the framework used to obtain the hard-thermal- loop self-energy in a system 
with an anisotropic momentum space distribution. In Sec. lIIII we present a tensor decomposition of 
the self-energy and dielectric tensors. In Sec. lIVI we work out the details of the tensor decomposition 
and give expressions for the self-energy "structure functions." In Sec. we discuss the static limit 
of the various self-energy structure functions. In Sec. IVII we use the tensor decomposition of the 
dielectric tensor to determine dispersion relations for all stable and unstable modes. In Sec. lVlll we 
present analytic expressions for the self-energy structure functions in the small-anisotropy limit. 
Finally, in Sec. IVII II we present conclusions and outlook for application of the results found here. 
We provide a summary of our notational conventions and expressions for the various self-energy 
structure functions in two appendices. 



II. HARD-THERMAL-LOOP SELF-ENERGY 



We begin by repeating some of the steps necessary to derive the hard-thermal-loop resummed 
gluon self-energy within semi-classical transport theory |^. Within this approach partons are 
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described by their phase-space densities and their time evolution is given by Vlasov-type transport 
equations 0, 0|- In this paper we will concentrate on the physics at the soft scale, k ~ gT <C T, 
which is the first scale at which collective motion appears. At this scale the magnitude of the field 
fluctuations is ^4 ~ ^/gT and derivatives are of the scale dx ~ gT. With this power-counting a 
systematic truncation of the terms contributing to the transport equations for soft momenta can 
be realized. 

At leading order in the coupling constant the color current, J^, induced by a soft gauge field, 
A^^, with four-momentum K = (w, k) can be obtained by performing a covariant gradient expansion 
of the quark and gluon Wigner functions in mean-field approximation. The result is 

&X)=9j ^V'^i2N,SN-ip,X)+Nfi6nl{p,X)-Snl{p,X)) , (1) 

where V^^ = (l,k/a;) is the gauge field four-velocity, 5N"-{p,X) is the fluctuating part of the 
gluon density, and 5n\{p, X) and 5n'^{p, X) are the fluctuating parts of the quark and anti-quark 
densities, respectively. Note that SN"" transforms as a vector in the adjoint representation {6N = 
5N°'T°') and dnfi^ transforms as a vector in the fundamental representation {5n± = 6n'^t°'). 
The quark and gluon density matrices above satisfy the following transport equations 

[V-Dx,Sn±{p,X)] = TgV^F>'''{p,X)d,n±{p) , (2) 
[V-Dx,6N{p,X)] = -gV^F^'''{p,X)d,N{p) , (3) 

where Dx = dx + igA(X) is the covariant derivative. 

Solving the transport equations and (jSJ for the fluctuations 5N and Suj- gives the induced 
current via Eq. 

J-UX) =9' J ^V^V%^f{p) J dr U{X, X - Vr)F^p{X - Vr)U{X - Vr, X) , (4) 
where U{X, Y) is a gauge parallel transporter defined by the path-ordered integral 

rY 



U{X, Y) = V exp 



Jx 



(5) 



Fai3 = daAp — dpAa — ig[An,Ay] is the gluon field strength tensor, and 

/(p) = 2iVeiV(p) + Nf{n+{p) + n_(p)) . (6) 

Neglecting terms of subleading order in g (implying U ^ 1 and F^p daAp — dpAa) and per- 
forming a Fourier transform of the induced current to momentum space we obtain 

W = i j ^V%>f(v) - J^^) A-'(K) , (7) 

where e is a small parameter that has to be sent to zero in the end. 

From this expression of the induced current the self-energy is obtained via 

which gives 
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This tensor is symmetric, Il'^''{K) = I['^'^{K), and transverse, Kf^Ilf^'^{K) = 0. Note that the same 
result can be obtained using diagrammatic methods if one assumes that the distribution function 
is symmetric under p — > — p 0|. 

In the hnear approximation the equation of motion for the gauge fields can be obtained by 
expressing the induced current in terms of the self-energy 

J.UK) = Il^''iK)A'^iK), (10) 

and plugging this into Maxwell's equation 

- zK^F^^'^iK) = Jr„d W + Jext W , (11) 

to obtain 

[K^g^''' - K^'K" + Yl^''' {K)]A^{K) = -J^^t{K) , (12) 

where is an external current. Using the gauge invariance of the self-energy we can write this 
in terms of a physical electric field by specifying a particular gauge. In the temporal axial gauge 
defined by = we obtain 

[(/fc2 - uj^)S^i - k'k^ + n'^{K)]E^{K) = {A-\K)y^E^{K) = iuj JL(i^) . (13) 

Inverting the propagator allows us to determine the response of the system to the external source 

E\K) = iu;A^HK)JUK). (14) 

The dispersion relations for the collective modes can be obtained by finding the poles in the 
propagator A^^{K). 



III. TENSOR DECOMPOSITION 



In this section we develop a tensor basis for an anisotropic system in which there is only one 
preferred direction. As mentioned above the self-energy is symmetric and transverse. As a result 
not all components of 11'^'^ are independent and we can restrict our considerations to the spatial 
part of U^'^, denoted II'-'. We therefore need to construct a basis for a symmetric 3-tensor that - 
apart from the momentum k^ - also depends on a fixed anisotropy three- vector n*, with = 1. 
Following Ref. P| we first define the projection operator 

j^ij ^ §ij _ k^k^/e, (15) 

and use it to construct n* = A^^n^ which obeys n • A; = 0. With this we can construct the remaining 
three tensors 

B'^ = k'y/k^ (16) 

C'^ = n^ni/h^ (17) 

D'^ = k'n^ + k^h\ (18) 

Any symmetric 3-tensor T can now be decomposed into the basis spanned by the four tensors 
A,B,C, and D 

T = aA + 6B + cC + (iD . (19) 

Furthermore, the inverse of any such tensor is then given as 

1 _ 1 (a + c)B-a-i(6c-nW)C-dD 
^ ^+ 6(a + c)-nW • ^^^> 
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IV. SELF-ENERGY STRUCTURE FUNCTIONS 



The spacelike components of the self-energy tensor can be written as 

At this point the distribution function /(p) is completely arbitrary. In order to proceed we need 
to specify a form for the distribution function. In what follows we will assume that /(p) can be 
obtained from an arbitrary isotropic distribution function by the rescaling of only one direction in 
momentum space. In practice this means that, given any isotropic distribution function /iso(p^), 
we can construct an anisotropic version by changing the argument 

/(p) = /iso(p' + ^(p-fi)') , (22) 

where n is the direction of the anisotropy and ^ > — 1 is a adjustable anisotropy parameter. Note 
that ^ > corresponds to a contraction of the distribution in the n direction whereas — 1 < ^ < 
corresponds to a stretching of the distribution in the n direction. This assumption allows us to 
simplify (|21l) by performing a change of variables to p 

p'=p'(l + e(v-n)2) . (23) 

After making this change of variables it is possible to integrate out the |p|-dependence giving 

WHK) ^ r.1 f '^.'l±^ { + . (24) 



47r (1 + C(v.n)2)2 \^ K-V + ie^ 
where 

We can then decompose the self-energy into four structure functions 

n*^' = aA'^ + pB'^ + -fC'^ + SD'^ , (26) 
which are determined by taking the following contractions: 

fflV^h^ = n^(a + 7) , 
Trrf^' = 2a + /3 + 7. (27) 

In Appendix^ we collect the resulting integral expressions for the structure functions. All four 
structure functions depend on niD, io, k, ^, and k • n = cos0„. In the limit ^ — > the structure 
functions a and P reduce to the isotropic hard-thermal-loop self-energies and 7 and 6 vanish 

a{K,0) = Ut{K) , 

P{K,0) = ^Ul{K) , 
l{K,0) = 0, 

6iK,0) = 0, (28) 
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with 



Ut{K) 



UJ 



2 



2k ^ UJ 



2ujk 
UJ + k 



k"^ UJ + k 



UJ 



k 



(29) 



(30) 



For finite ^ the analytic structure of the structure functions is the same as in the isotropic case. 
There is a cut in the complex uj plane which we can chose to run along the real uj axis from 
—k<uj<k. For real- valued uj the structure functions are complex for all a; < A; and real for 
UJ > k. For imaginary- valued uj all four structure functions are real- valued. In Fig. ^ we plot the 
structure function a for real and imaginary values of ^ = {0, 1, 10}, and 6n = 7r/4. 

With these structure functions in hand we can construct the propagator A^^(K) using the 
expressions from the previous section. Writing A.~^{K) in terms of our tensor basis 

- + a) A + (/3 - uj^)B + -fC + 6U 



and applying the inversion formula H20|) we obtain an expression for the propagator 
A{K) = AaA + {k'^ - uj"^ + a + 7)AgB + [(/3 - uj'^)Ag - Aa]C - 5AgD 

with 



k' 



UJ + a 

,2 



Note that we can reorganize A{K) and write it as 

A{K) = Aa [A - C] + Ag [{k^ -uj^ + a + ^)B + {(3 



uj^)C 



5D] . 



(31) 

(32) 

(33) 
(34) 

(35) 



V. STATIC LIMIT 

In order to see how the momentum-space anisotropy in the distribution functions affects the 
response to static electric and magnetic fluctuations we examine the limit u; — > of the propagators 
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FIG. 2: Angular dependence of m^, ml, and ml at fixed (a) ^ = 10 and (b) ^ = —0.9. 



(jSnj and H34|) . Approaching along the real co axis we find that to leading order a ~ 7 ~ 0{uj^), 
13 ~ 0{uj'^), and 5 ~ 0(ia;).^ We can therefore define four mass scales 



m„ 



m^ 



ml 



lim a , 

lim 2 ' 

lim 7 , 

lim Im(5 



Writing the static limit of the propagators H33() and (|34() in terms of these masses gives 



a; 



UJ 



2 ^ 



A;2 



(A;^ +m^ +m^)(A;2 +m|) 



4 

m^ 



A(^^ can be factorized into 



UJ 



where 
with 



2m| = M^± JM^ - 4(m2 (m2 + rnl) 



m 



m. 



2 _ 



2 , 2 , 2 

m^, + m^ + m. 



(36) 

(37) 
(38) 

(39) 
(40) 
(41) 
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m| 



m 



73 ' ■'-7 

2—0 and m^ = m|). For finite it is possible 



In the isotropic limit, ^ 0, m^ 
to evaluate all four masses defined above analytically. The results for nia and m^ are listed in 
Appendix^ In Fig. I^Jwe plot the angular dependence of m^ , m^, and ml at fixed = 10 and 
^ = —0.9. In the case ^ > (Fig. [2^) we see that for small 9n the scale m^ > m^j and for On 
near 7r/2, m^ < m^. For small 0^ the scales m^ and ml are negative. The fact that rua and m_ 
are non- vanishing is in agreement with the findings of Cooper et al. ^3]; however, they neglected 



^ Identical results can be obtained by coming in along the imaginary axis with a suitable redefinition of rrig. 
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to consider the fact that these masses might be negative and would therefore not correspond to 
screening of the magnetic interaction. The fact that these quantities are negative indicates that 
for ^ > the system possesses an instability to transverse and "mixed" external perturbations 
associated with and m^, respectively. The transverse instability is present for any 0„ 7^ 7r/2 
while the mixed instability is only present for On < ^^"^'^ with 0'^^^'^'^ depending on the value of 
^. In the case ^ < (Fig. [2h) we see that for small On the scale < rnjj and for On near vr/2, 
> m^). For On ^ 7r/4 the scale is negative again signaling the presence of an instability in 
the system. In the next section we will discuss these instabilities in more detail. 



VI. COLLECTIVE MODES 

A similar factorization of A^^ can be achieved in the non-static case allowing us to determine 
the dispersion relations for all of the collective modes in the system. 



A. Stable Modes 

First, let's consider the stable collective modes which have poles at real-valued uj > k. In this 
case we factorize A^^ as 

A^'^ = {io^ - nl){u;^ - nl) , (42) 

where 

2nl = n'^± ^n^- 4((a + 7 + P)/? - k'^n^d'^) , (43) 

and 

J^2 = ct + /5 + 7 + A:2. (44) 

Note that the quantity under the square root in H43|) can be written as (a — /3 + 7 + /c^)^ + Ak'^'n?'5'^ 
which is always positive for real uo > k. Therefore there are at most two stable modes coming from 

The remaining stable collective mode comes from the zero of A^^. The dispersion relations for 
all of the collective modes can be determined by finding the solutions to 

= 0|(a;±) , (45) 
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FIG. 4: Ta{k) and r_(A;) as a function of /c with (a) ^ = 10 and 9n = vr/8 and (b) ^ = —0.9 and 
On = vr/2. 



= A:^ + a(u;Q,) . 



(46) 



In the isotropic hmit (|28|1 cJq, = = liJj- and w„ = cjj^. For finite there are three stable 
quasiparticle modes with dispersion relations which depend on the angle of propagation with respect 
to the anisotropy vector, In Fig. |31 we plot the dispersion relations for all three modes for 
mo = \/3, ^ = 10, and On = {0, 7r/4, 7r/2}. 

B. Unstable Modes 

For non-zero ^ the propagator also has poles along the imaginary uj axis.^ The dispersion 
relation for these modes can be determined by taking w — > iF with F real-valued and solving for 
F(A;). In this case we factorize the inverse propagator as 

- (F^ + 17^)(f2 + 1^2) ^ 



A 



G 



(47) 



where i7± on the right hand side are evaluated at w = iF. However, in contrast to the stable modes 
there is at most one solution in this case since numerically we find that Vl\ > Q for all F > 0. 

For ^ > there is also an unstable mode present in so that in this case there are two 
unstable modes in the system which can be found by solving 



-Jl_(iF_) , 
-k^ — aiiVr 



(48) 
(49) 



Note that in both cases there are two solutions corresponding to modes with positive and negative 
growth rates. One of these corresponds to an exponentially growing solution and the other an 
exponentially decaying one. In Fig. 0^ we plot Ta{k) and F_(A;) with ^ = 10 and On = it/8. For 
^ < there is no longer an unstable mode coming from A^ and there is therefore only one unstable 
mode coming from F_. In Fig. we plot F_(/c) with ^ = —0.9 and On = 7r/2. 



We have checked for poles at complex lu numerically but found none. 
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VII. SMALL ^ EXPANSION 



In the small-^ limit it is possible to obtain analytic expressions for all of the structure functions 
order- by-order in ^. To linear order in ^ 

— (3 5 cos 26'„)m|, - -{1 + cos26'„)m|, 



-2/3 = Ul{z)+^ 



-^Uriz) ((1 3 cos 26i„) - z^ (3 + 5 cos 

^2 



29n 



1 Z^ 

-{l + 'ims29n)m]. + IlL{z) coslOn (1 + 3 COS 26*^ 

6 \ 2 
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= -^{4z^ml + 3UTiz)il-4z^))cos9n 



where z = u/k. 



A. Static Limit 



(50) 



Using the linear expansions and the fact that in the static limit 11^ ~'"T'D ^'^'^ 
-muj/{Ak) we can write for the masses ()36() 



VTLf. 



-|(l + cos2e„) , 
b 

l + |(3cos20„- 1) , 
o 

- sm dn , 
vr 

-^-sin6'„cos6'„ , 



(51) 



where m = m /mj-). Using these we can obtain small-^ expressions for m± defined in (|4U() 



m_ 



f 

1 + |(3cos26l„ - 1) , 



6 

cos 29r, ■ 



(52) 



B. Collective Modes 



Since 5 is O(^) it can be ignored in the expansion of (|43() so that to linear order in ^ the collective 
modes satisfy 



-uj^ + a = 



(A:2-a;2 + a + 7)(/?-c^2) = 0, 



(53) 



where a, /3, and 7 are given by (|5U|) . Note again that there is only one unstable mode coming from 
A^^ since /3(ir) > for all T > 0. 
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VIII. CONCLUSIONS 



In this paper we have derived a tensor basis for the gluon self-energy in a high-temperature 
quark gluon plasma with an anisotropic momentum-space distribution. We then restricted the 
distribution function by requiring that it could be obtained from an isotropic distribution function 
by the rescaling of one direction specified by an anisotropy vector, fi, and strength, ^. Positive 
values of ^ correspond to a contraction of the isotropic distribution function along n while negative 
values of ^ correspond to a stretching along n. Within this framework we could derive analytic 
forms for all of the structure functions associated with the tensor basis. Using these expressions 
we were then able to identify and determine the dispersion relations for the collective modes for 
both positive and negative ^. We found that for ^ > there were at most three stable and two 
unstable modes with dispersion relations which depended on the angle between the wave vector, 
k, and the anisotropy vector. For ^ < we found that there were also three stable modes but only 
one unstable mode. Additionally, we obtained analytic expressions for the structure functions in 
the limit of small ^. These results should provide a reference point for the systematic study of the 
isotropization of a relativistic plasma. 

The study of Mrowczyhski and Randrup suggests that during heavy-ion collisions the rate of 
isotropization via collective modes is comparable with collisions and therefore cannot be ignored 
In this paper we have made no attempt to discuss the phenomenological rate for instability 
growth because there are a number of questions which would need to be addressed prior to making 
any definitive statements about the role of instabilities in plasma evolution and their expected 
contribution to observables. This is because we have only derived the self-energy in a linear 
expansion in the fluctuations and to leading-order in the coupling constant. Assuming that there 
is truly exponential growth of the fields in the direction of the anisotropy this means that the 
linear approximation will break down very quickly. In practice the non-linear terms in the transport 
equations will become important and regulate the growth of the modes which have become unstable. 

Within electrodynamics the coupling constant is small and it is possible to experimentally 
study the Weibel instability However, with QCD the story is dramatically different since 

for experimentally realizable situations the coupling constant is large and the non-linear effects 
due to gluon self-interaction become important much sooner than any non-linear effects would for 
QED. Nevertheless, this does not diminish from the fact that these unstable modes exist and will 
therefore have a role to play in plasma evolution. In order to assess this role, however, detailed 
studies of the time evolution of anisotropic quark-gluon plasmas will need to be performed. 
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APPENDIX A: NOTATION AND CONVENTIONS 

We summarize here the notation and conventions which we use in the main body of the text. 
• Natural units: h = c = 1 
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Metric: g^''' = g^y = diag(l, -1, -1, -1) 

4- vectors: Indicated by Greek indices, e.g. K^^ = (a',k). 

3- vectors: Indicated by lowercase Latin characters. Upper Latin indices like A; use a 
Euclidean 3-metric, e.g. k = fc*, = k^. 



Fourier transform: 



3 

f{X) 



AK X 



f{K) 



APPENDIX B: ANALYTIC EXPRESSIONS FOR STRUCTURE FUNCTIONS 



In this appendix we collect the integral and analytic expressions for the structure functions a, 
/3, 7, and 5. Choosing n = z and k to lie in the x—z plane {kx/kz = tan On) we have v • n = cos and 
V • k = kxCos(f)sm.9 + kz cos 9. Using this parameterization the (j) integration in all four structure 
functions defined by the contractions in Eq. ()27() can be performed analytically. 



f d{cosO) u + ^kz cos 6 

k'^n^ 



2 (1 +^(cos 0)2)2 

+k'^{s^ - (cos 



00 — kz cos i 



U)kz ,2 



) )R{uj — kz cos 9, kx sin ( 



(B.l) 



mjjuj^ f d{cos 9) 



1 



k^ J 2 (l+e(cos 0)2)2 
^(cos^) 1 



[1 - (lj + (,k, cos 9)R{uj- kz cos 9, kx sin 9)] , (B.2) 



2 fc2(l +^cos2 0)2 
{uj + ^kz cos 9)k^ 



k^ 



up" + ik^ cos^ - - ^k^ cos^ 9) 



1.2 



2 (cos (9 



u)k 



z \2 



A;2 



I R{uj — kz cos 0, kx sin 9) 



,(B.3) 



r/r^ ^^ rn'rjUJ f d(cos 9) ui + S,kzCos9 / ,,2 ^ , n r^/ , /, , ■ /t^ 

S{K,0 = \ n ,1 ^—^ ik, + {k^cose-ujkz)R{uJ-kzCOs9,kxSme)) , B.4 

k^w^J 1 (l+4cos^yj^ ^ ' 



where s2 = (1 — u? Ik^){k\lk^) and 

S-2-K 

R{a,b) 



lo 2tt a - b cos (j) + ie \/a + b + ie\J a - b + ie 
When a and b are real-valued R can be simplified to 

sgn(a)e(a2 - b2) ie(62 - a^) 



R{a,b) 



- 62 



(B.5) 



(B.6) 



with 0(x) being the usual step- function. Note that the remaining integration over 9 can also be 
done analytically but the results are rather unwieldy so we do not list them here. 
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STATIC LIMIT 



In the hmit a; —> it is possible to obtain analytic expressions for all four structure functions. 
The results for rua and mp defined in Eq. are 



m 



,2 



m 



2 



A;^ arctan ^ - arctan ^ = I , (B.7) 



+ (1 + ^)arctan^)(k2 + ^k^) + ^k, f k,^ + Ji^arctan^= 



with similar results for and m|. 
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